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Abstract 

We study effects of nonlocality in the nuclear force on the G-matrix elements for finite 
nuclei. Nuclear G-matrix elements for 16 O are calculated in the harmonic oscillator basis 
from a nonlocal potential which models quark exchange effects between two nucleons. We 
employ a simple form of potential that gives the same phase shifts as a realistic local 
nucleon potential. The G-matrix elements calculated from the nonlocal potential show 
moderate increase in repulsion from those derived from the local potential. 



1 Introduction 

Several models of the nucleon-nucleon interaction have been presented and applied to nuclear 
physics problems. The G-matrix interaction plays a central role in calculations of nuclear 
structure from the nucleon-nucleon (NN) interaction, as it is the basic two-body interaction in 
many body systems. 

It is well known that the nucleon-nucleon interaction is described by the quark cluster model 
(QCM). This model gives a good description of the nucleon-nucleon scattering |], |2j. The main 
feature of the interaction is nonlocality at short distances. This short range potential has a 
nonlocal Gaussian soft core of the form, 

r — r' 

V(r', r) oc exp[-(— — ) 2 ]. 

This type of nonlocality is expected from the quark antisymmetrization between the baryons. 
The Gaussian form reflects the quark wave function inside the baryon and the nonlocality range 
parameter b is determined by the size of the baryon. The nonlocal exchange force can reproduce 
the short-range repulsion of the nuclear force at low energy. The repulsive "core" in QCM is 
a soft core with energy dependence in its equivalent local form. It allows the NN relative 
wave function to go inside and thus gives a milder form factor of the deuteron (or nuclei) at 
large momenta || . Because the quark exchange nonlocality comes only at short distances, the 
central part of the potential is most affected. In fact, the NN repulsion is roughly spin-isospin 
independent. 

In this paper, we would like to study significance of the nonlocality of quark model origin 
in nuclear structure calculations. Since we can not see signals of such nonlocality in on-shell 
properties, it would be interesting if we could see them in nuclear phenomena. To do so, we 
calculate G-matrix elements from nonlocal potentials. In conventional nuclear physics, the G- 
matrix interaction is used as the basic two-body interaction for the nuclear many-body problem. 
Hence, if there would be signals of quark-originated nonlocality in nuclear phenomena, they 
must exist in the G-matrix interaction as well. Previously, effects of symmetry restoration of 
QCD in G-matrix elements were examined Q. Thus in the present study we investigate another 
possibility of QCD oriented effects in nuclear physics. 

This paper is organized as follows. In section 2, a method to calculate G-matrix elements 
for finite nuclei is given briefly. In section 3, we introduce the Tamagaki G3RS (Gaussian 
soft core potential with three ranges) potential [0] as a local potential. Then we construct a 
nonlocal Gaussian soft core potentials, which has the same intermediate and long range parts 
as the original local potential. The short-range part is replaced by a nonlocal potential with 
free parameters which are determined so as to reproduce the same scattering phase shifts as 
the local potential. In section 4, numerical results for the G-matrix elements are given both for 
the local and nonlocal potentials. A summary is given in section 5. 

2 Nuclear G-matrix for Finite Nuclei 

The G-matrix is the most fundamental two body interaction which incorporates the simplest 
many body effect, i.e., the Pauli blocking. It is obtained by solving the Bethe-Goldstone 
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equation ||: 

G(E)=v + v ¥ ® w G(E). (1) 

Here v is a nucleon-nucleon interaction in free space, Q the Pauli-blocking operator, and H 
a single particle hamiltonian. We compute G-matrix elements for finite nuclei, specifically for 
16 0, following the method developed in previous works []7|, |8], |9|. Hq is chosen to be the harmonic 
oscillator hamiltonian. The two-body problem is then reduced to a one-body problem of the 
relative coordinate by employing the Eden and Emery's approximation |TD[ for the Q operator. 
Explicitly, after reducing the Pauli-blocking operator as a function of the relative coordinate, 
it is assumed to be 

O(n') = \ 1 for n ' > n + 2 (2) 
I otherwise 

where n is the node quantum number for the relative motion of the initial two nucleons and 



n' that for intermediate states. The Barret-Hewitt-MacCarthy's method [IT] is then used for 
matrix inversion. Additional parameter is introduced for the gap energy, which is chosen to be 
40 MeV to achieve a good agreement between the calculated G-matrix elements and empirical 
ones ||12|| . Finally, the starting energy is averaged over the Fermi surface by assuming that the 
relevant two nucleons are there. 

G-matrix elements are obtained first in partial wave channels denoted by 1 S'o, 1 Pi, etc., 
which are then transformed into those in interaction channels denoted by SE (Singlet Even), 
TE (Triplet Even), etc. Matrix elements in SE and SO (Singlet Odd) channels are just those in 
1 Sq and X P\ channels, respectively. The TE and TNE (Tensor Even) components are obtained 
from the coupled 3 Si — 3 D\ channels: 

G(TE) = (S | G( 3 Si - 3 £>i) | S), (3) 
G(TNE) = (S \G( 3 S 1 - 3 D ± ) \ D). (4) 

The TO, TNO, LSO (LS Odd) and LSE (LS Even) components are obtained from the following 
relations 0: 

G(TO) = G{ 3 P ) + 2G{LSO) +AG{TNO), (5) 
G(TNO) = -A[ 2G (3p o) _ 3G( 3p i) + G( 3p 2)]) (6) 

G(LSO) = -^[2G( 3 P ) + 3G( 3 P 1 )-5G( 3 P 2 )], (7) 
G(LSE) = -L[9G( 3 D 1 )+5G( 3 D 2 )-UG( 3 D 3 )], (8) 

In the previous studies [0, ^ |9|, the Reid soft-core potential, Paris potential, and other 
phenomenological nuclear forces have been applied. We here choose the Tamagaki G3RS po- 
tential m, which is a three range gaussian with a short range soft core. We replace the shortest- 
range part of the original potential by a nonlocal one, and apply it to the G-matrix calculation. 
The resulting G-matrix elements are compared with each other. 
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3 Gaussian Nonlocal Potential 



We use the following Tamagaki potential as a standard local potential, 

V L (r) = V c {r) + S 12 V T (r) + (L ■ S)V LS (r) + W 12 V w (r) + L 2 V LL (r), (9) 

with 

S 12 = 3(<7i • r)(cr 2 • r)/r 2 - o"i ■ a 2 

W l2 = 1/2{(<7i • L)(<7 2 - L) + (cr 2 • L)(<7x ■ L)} — (oi • cr 2 )L 2 /3 
= (L-S) 2 -{5 LJ + (a 1 - ( x 2 )/3}L 2 

where L and J stand for the orbital and total angular momenta, respectively. The radial 
functions V n {r) (n = C, LS, . . . ) are given by 

V n (r) = j2Vmexp[-(r/ Vni )% (10) 

i=l 

where the suffices i = 1,2,3 refer to the short range, intermediate range and long range part, 
respectively. The potential parameters, V n i and r] ni , are determined so as to reproduce the 
scattering phase shifts as given in Table [l] || . 

The nonlocal potential is constructed by replacing the short range part (i = 1) of Eq. (JT0|) . 
We consider nonlocality only for the central part of the form: 

V(v',r) = V NL (* 1/2 b)- 3 ex P [-(^) 2 - (^) 2 ] (11) 

Similar nonlocal terms might appear in noncentral forces, but for the reason explained in 
introduction, we consider the nonlocality only in the central potential. In Eq. (|TT|), the factor 
(n 1 / 2 ^^ 3 exp[— ( L j L ) 2 } represents the nonlocality of the nuclear force at short distances, with 
a new parameter b as the range of nonlocality. When & — 0, this part reduces to the delta 
function, 5(r — r'), and then the nonlocal potential becomes the local one. In the quark cluster 
model, b is proportional to the size of the quark wave function in the nucleon. A typical value 
of b would be 0.5 fm, which we use in the following calculations. 

In order to study effects of nonlocality exclusively, we require the nonlocal potential to 
reproduce the same scattering phase shifts as the corresponding local potential does. We call 
such a nonlocal potential as local-equivalent (LE) nonlocal potential. For this, we adjust the 
potential parameters, Vjvl an d i]nl- We allow them to depend on the angular momentum L. 
Phase shifts in various partial waves at scattering energies from E cm = to 300 MeV are fitted. 
As shown in Figs. 1-2, both the local and LE nonlocal potentials (LENL) reproduce empirical 
phase shifts well up to the scattering energy ~ 300 MeV. In particular, the difference between 
the local and LENL are negligibly small. In Figs. 1 and 2, we also present the phase shifts 
calculated from the nonlocal potentials (NL) with the same strength, Vnl — Va and the same 
range t/jvl = Vci as the original Tamagaki potential. They show that the nonlocality yields a 
softer repulsion at large energy. Although we have shown here the phase shifts only for the two 
channels, we find qualitatively the same feature for other S and P channels. 
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Table 1: Parameters of the Tamagaki G3RS potential 
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Table 2: Parameters of the nonlocal potential with b = 0.5 fm 

Summarizing, our LE nonlocal potential reads 

V NL (r\v) = ^ Vi (vr 1 / 2 6)- 3 exp[-(^) 2 -(^) 2 ] 

+ Va exp[-(r/ Vci ) 2 ] 6(t - r') (12) 

with the parameters given in Table [l] for the local part and in Table || for the nonlocal part 
(b = 0.5 fm). 

4 Numerical Results 

The G- matrix elements in the harmonic-oscillator basis are presented in Tables [3|-[| for the local 
(Tamagaki) and LE nonlocal potentials with b = 0.5 fm as functions of node quantum numbers 
n and n'. These are calculated using an oscillator parameter huo = 14 MeV for 16 and with 
the gap energy of 40 MeV, which are the same as those used in Refs. J7| ||, |[. 

First we compare G-matrix elements calculated from various nuclear forces. For this pur- 
pose, we show in Table [5] the G-matrix elements derived from the Paris potential and those of 
the present calculation for the SE channel. It turns out that both matrix elements are very 
similar. The difference is significant only for less important off-diagonal matrix elements, which 
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Table 3: G-matrix elements from the local potential for 16 
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Table 4: G-matrix elememnts from LE nonlocal potential (b=0.5fm) for 16 
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Table 5: Comparison of diagonal G-matrix elements in the SE channel derived from three 
nuclear forces. 

is, however, typically 20 % level. This is interesting because the Tamagaki potential has the 
Gaussian tail and its functional form differs from the others. This fact suggests that G-matrix 
elements are strongly constrained by on-shell properties (phase shifts). We have confirmed 
this by calculating G-matrix elements using NN potentials which do not necessarily reproduce 
empirical phase shifts. 

Now we turn to see the effects of nonlocality in the G-matrix elements. From Tables |3|-[|, 
one sees that for non-central channels of tensor and LS, the results of the local and of LE 
nonlocal potentials are almost identical. This seems natural because in the present calculation 
nonlocality is introduced only in the central channel. 

Effects of nonlocality in the G-matrix elements are best seen in the S-wave channels of SE 
and TE. They are repulsive, though the effects are typically as small as a few % or even less. 
Therefore, one may conclude that the present nonlocality from the quark exchange effects are 
negligibly small as compared with other many-body effects of nuclear physics which is typically 
about 10 % or even more [|l!|. The reason that the present nonlocality acts as repulsive 
is understood in the following way. As we have discussed in section 3, when nonlocality is 
introduced for a repulsive component, it effectively reduces the repulsion in the NN potential. 
To get the local-equivalent phase shifts, one has to modify the potential parameters by reducing 
the range and increasing the strength of repulsion as summarized in Table [2| In the calculation 
of the G-matrix elements, the Pauli-blocking operator Q forbids the transition to lower levels 
which are already occupied. This effectively reduces the chance for the two nucleons to feel 
attraction in the nuclear force which comes mainly from the transition to the lower states. 
Hence, in the G-matrix elements, repulsive components are more enhanced. 

The nonlocal effect would become important for heavier nuclei where more levels are Pauli 
blocked. We have performed a calculation of G-matrix elements in the zirconium region. The 
calculation was done in the same way as for 16 0, but with tkv = 8.8 MeV and with the following 
Pauli-blocking operator: 

n(n') = { 1 for ™ />n + 4 ' (is) 
[ otherwise . 

Here, as the mass number is increased, more states are Pauli-blocked, which is implemented in 
the inequality n' > n + 4 in place of n! > n + 2 for 16 in (g). Results are summarized in Table 
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Table 6: G-matrix elements in the SE channel for zirconium 

|5] for the SE channel, where we see slightly more nonlocal effects. The difference is, however, 
once again very small. 

5 Summary 

The G-matrix elements for 16 have been calculated in the harmonic oscillator basis from 
a G3RS Tamagaki potential (local potential) and a nonlocal potential which has a simple 
Gaussian nonlocality at a short range but produces the same scattering phase shifts as the 
local potential does. The nonlocal Gaussian soft core is typical of quark cluster model (QCM) 
which gives a good description of the nucleon-nucleon scattering. 

The G-matrix elements thus derived from the nonlocal potential show an effective repulsion 
in comparison with those derived from the local potential. The difference appears in the SE 
and TE channels mostly. The effects are, however, generally very small for all nuclei. It is 
not, therefore, likely that quark-originated nonlocal effects are detected in conventional nuclear 
phenomena. 
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